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The disorder induced metal-insulator transition is investigated in a three-dimensional simple cubic 
lattice and compared for the presence and absence of time-reversal and spin-rotational symmetry, i.e. 
in the three conventional symmetry classes. Large scale numerical simulations have been performed 
on systems with linear sizes up to L = 100 in order to obtain eigenstates at the band center, 

E = 0. The multifractal dimensions, exponents D q and a q , have been determined in the range of 
— 1 < q < 2. The finite-size scaling of the generalized multifractal exponents provide the critical 
exponents for the different symmetry classes in accordance with values known from the literature 
based on high precision transfer matrix techniques. The multifractal exponents of the different 
symmetry classes provide further characterization of the Anderson transition, which was missing 
from the literature so far. 
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I. INTRODUCTION 

The metal-insulator transition (MIT), and disordered 
systems have been at the forefront of condensed matter 
research since the middle of the last centurj®, and yet 
this topic still has several open questions and is still ac¬ 
tively investigated. In the last few years experimental 
evidence has been obtained about this topic; in particu¬ 
lar, reporting Anderson localization of ultrasound in dis¬ 
ordered elastic network ^ 2 light in disordered photonic 
lattices in the transverse direction®, or in an ultracold 
atomic system in a disordered laser trapP. Richardclla 
et alW examined the MIT in a dilute magnetic semi¬ 
conductor Gai-^-Mn^As, which is a strongly interacting 
and disordered system. They found a clear phase tran¬ 
sition together with multifractal fluctuations of the focal 
density of states (LDOS) at the Fermi energy, showing, 
that multifractality is a robust and important property 
of disordered systems. Multifractal properties consistent 
with the theory of Anderson localization are also found 
in the ultrasound system 3 . On the theoretical side, we 
know that disorder plays a crucial role in integer quantum 
Hall effect^, and recently it was shown that an enhanced 
correlation of multifractal wave-function densities in dis¬ 
ordered systems can increase the superconducting criti¬ 
cal temperature® or the multifractal fluctuations of the 
LDOS close to criticality may lead to a new phase due 
to the presence of local Kondo effects induced by focal 
pseudogaps at the Fermi energjP. Moreover, Anderson 
localization has also been reported in the spectrum of the 
Dirac operator within the lattice model of QCD at high 
temperatures using spectral statistics 1 ^, and multifractal 
analysis seems to corroborate it, as well 11 . 

These models show an increased interest in under¬ 
standing the nature of the Anderson transition in the 
presence of various global symmetries. A comprehensive 
review of the current understanding is given in Ref. 1121 
These symmetry classes have been introduced first to de¬ 
scribe random matrix ensembles, but the naming con¬ 


ventions are the same in the field of disordered systems. 
The classification considers two global symmetries: time- 
reversal and spin-rotational symmetry. As it turns out, 
beside these symmetries there are three further symme¬ 
try classes according to the presence of chiral symme¬ 
try, and in addition there are four Bogoliubov-de Gennes 
classes also, corresponding to particle-hole symmetry 12 
prominent in hybrid (superconductor-normal) systems. 
The effect of symmetry classes at the Anderson transi¬ 
tion has already been investigated earlier 1 ® using spec¬ 
tral statistics, but there is much less work based on the 
multifractal analysis of the eigenstates, and multifractal 
exponents are known numerically only for the orthogonal 

clasJi®. 

Our goal in this article is to fill in this gap and apply 
multifractal finite-size scaling (MFSS), developed orig¬ 
inally by Rodriguez, Vasquez, Romer and SlevirP®, to 
the Anderson models in the three conventional Wigner- 
Dyson (WD) classes. The organization of the article is 
the following. In Sec. |TT] we define the model and de¬ 
scribe its numerical representation. In Sec. |III| we briefly 
describe the finite-size scaling analysis of the general¬ 
ized multifractal exponents of the critical eigenstates, in 
Sec. [IV] we give the results obtained for the three uni¬ 
versality classes and finally in Sec. [V]we summarize our 
results. 


II. MODELS AND NUMERICAL 
REPRESENTATION 

A. The model 

In this article we investigate Anderson models belong¬ 
ing to the three WD classes, without chiral and particle- 
hole symmetry. We investigate the case of diagonal disor¬ 
der and nearest-neighbor hopping, therefore the Hamil- 
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tonian reads as 


^ ^ ] e i c J CT c icr 

io 


^ 'j tijcra lC icr C jcr' h.C., 
ijaa' 


(i) 


where i, j and a, a' stand for site- and spin index, £j-s are 
random on-site energies, which are uniformly distributed 
over the interval [— , W acts as disorder strength. 

Using a uniform distribution is just a convention, other 
distributions of disorder, e.g. Gaussian, binary, etc. can 
be used as well. 

In the orthogonal class time-reversal and spin- 
rotational symmetry are preserved. In this case the 
Hamiltonian is invariant under orthogonal transforma¬ 
tions - hence the name therefore it is a real symmetric 
matrix. Since spin does not play a role, we consider a 
spinless Anderson model. In the numerical simulations 
the Hamiltonian is represented by an N x N real sym¬ 
metric matrix, where N = L 3 , and L is the linear system 
size in lattice spacing. The diagonal elements, are are 
uniformly distributed random numbers, the off-diagonal 
elements are zero, except if i and j are nearest neighbors: 


H?; = 


£i&U [-^, ^] , if i=j 
— 1 , if * and j are neighboring sites 
0 , otherwise 


( 2 ) 


The energy unit is fixed by setting the hopping elements 
to 1. To avoid surface effects, we use periodic bound¬ 
ary conditions. However, this case was investigated very 
carefully by Rodriguez et aZP^l, we consider this symme¬ 
try class to verify our numerical method, and to obtain 
a complete description of all the WD classes. 

In the unitary class time-reversal symmetry is broken, 
which can be realized physically by applying a magnetic 
field. It can be shown, that either spin rotational sym¬ 
metry is broken or not, the model will belong to the uni¬ 
tary clasJI-i The Hamiltonian is invariant under unitary 
transformations therefore it is a complex hermitian ma¬ 
trix. We discuss the case when spin-rotational symmetry 
is present, because this way we can use spinless fermions 
again, which keeps the matrix size N x N. However, one 
has to store about twice as much data compared with 
the orthogonal case, because here every off-diagonal ma¬ 
trix element is a complex number. Obviously finding an 
eigenvalue and an eigenvector takes more time, too. 

For the numerical simulations we followed Slevin and 
OhtsukPl. Let us consider a magnetic field pointing in 
the y direction with flux <h, measured in units of the flux 
quantum, h/e. Its effect can be represented by a unity 
phase factor, the Peierls substitution for the hopping el¬ 
ements of the Hamiltonian matrix. The upper triangular 


of the Hamiltonian reads as 

' ei eC/[-f,f],if* = j 

— 1 , if i and j are neighboring sites 
ij in the x or y direction 

_ e *27r$a^ ^ j and j are neighboring 
sites in the £ direction 
, 0 , otherwise 


( 3 ) 


Complex hermiticity sets the off-diagonal elements in the 
lower triangular part, j < i. Periodic boundary con¬ 
ditions and flux quantization force a restriction for the 
magnetic flux namely, that $ • L must be an integer. 
In the thermodynamic limit arbitrarily small magnetic 
field drives the system from the orthogonal to the uni¬ 
tary class. However, in a finite system the relationship 
between the system size, L, and the magnetic length, 
Lh = matters. In the case of weak magnetic field, 
L <C Lh, the system belongs to the orthogonal class, in 
the case of strong magnetic field, L 3> Lh , it belongs to 
the unitary class. Since we use system sizes that are mul¬ 
tiples of 10 lattice spacings, see Tab. |IJ we chose $ = |. 
This leads to Lh ~ 0.892 therefore this choice clearly 
fulfills the two conditions above. 

In the symplectic class time-reversal symmetry is 
present, and spin-rotational symmetry is broken, which 
describes a system with spin-orbit interaction. In this 
case the Hamiltonian is invariant under symplectic trans¬ 
formations therefore it is a quaternion hermitian matrix. 
For the numerical simulations we followed Asada, Slevin 
and Ohtsukf 16 . Since in this case we have to deal with 
the spin index also, the Hamiltonian is an 2 N x 2 N com¬ 
plex hermitian matrix. Diagonal elements corresponding 
to the 'ttli site and hopping elements between sites i and 
j are 2 x 2 matrices because of the spin indexes, having 
a form 

/ e . 0 \ ( e iaij cos P ij e ijij sin \ 

— I 0 £, ) ^ij = ( _ e -n« sinftj e -i<Xij cos p tj J > 

(4) 

where £* is an uniformly distributed random on-site en¬ 
ergy from the interval [— ^], a^, fiij and 7 ^ were 

chosen to form an SU(2)-invariant parametrization, lead¬ 
ing to the so-called SU(2) model: oy ? and 7 ^ are uniform 
random variables from the interval [0, 2-7t], and /? has a 
probability density function p(f3)d/3 = sin(2/3 )d/3 in the 
range [0, . The upper triangular of the Hamiltonian 

has the following form: 

Uu if 1=3 

H ?<i = tij , if i and j are neighboring sites (5) 
I 0 , otherwise 

The off-diagonal elements are defined following complex 
hermiticity. To store the Hamiltonian requires about 
eight times more memory compared to the orthogonal 
case, because here every off-diagonal element contains 
four complex numbers. Finding an eigenvalue is much 
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system size (L) 

number of samples 

20 

15000 

30 

15000 

40 

15000 

50 

15000 

60 

10000 

70 

7500 

80 

5000 

90 

4000 

100 

3500 


TABLE I: System sizes and number of samples for the simu¬ 
lation for each WD symmetry class. 


slower than for the unitary case, mainly because of the 
linear size of the matrix is twice as large. 


B. Numerical method 

MFSS deals with the eigenvectors of the Hamiltonian, 
which is a large sparse matrix. Recent high precision 
calculation^^ use Jacobi-Davidson iteration with incom¬ 
plete LU preconditioning, therefore we decided to use this 
combination. For preconditioning the ILUPACKpSl was 
used, for the JD iteration the PRIMMlPS package was 
used. Since the metal-insulator transition occurs at the 
band center 1 ^ {E = 0) at disorder W c ° ss 16.5 for the 
orthogonal, at W 1 / ss 18.3 for the unitary (depending 
on the strength of magnetic field), at W r f ss 20 for the 
symplectic class (for our parameters), most works study 
the vicinity of these points. To have the best compari¬ 
son, we analyzed this regime, therefore 20 disorder values 
were taken from the range 15 < W < 18 for the orthogo¬ 
nal class, 23 disorder values were taken from the interval 
17 < W < 20 for the unitary class, and 20 disorder val¬ 
ues were taken from the interval 19.4 < W < 20.5 for 
the symplectic class. System sizes were taken from the 
range L = 20.. 100, and the number of samples are listed 
in Tab. [TJ We considered only one wave-function per re¬ 
alization, the one with energy closest to zero in order 
to avoid correlations between wave-functions of the same 
system^. 


III. FINITE SIZE SCALING LAWS FOR 
GENERALIZED MULTIFRACTAL EXPONENTS 

In recent high precision calculationiSl the multifractal 
exponents (MFEs) of the eigenfunctions of the Hamilto¬ 
nian have been used to describe the Anderson metal- 
insulator transition. We use almost the same notation 
and methods as Ref. m but for better understanding 
here we introduce shortly the most important quantities 
and notations. The method has recently been success¬ 
fully extended for the investigation of the quantum per¬ 
colation transition in three dimension^ 17 ! 


Considering a d-dimensional cubic lattice with linear 
size L , one can divide this lattice into smaller boxes with 
linear size £. If 4/ is an eigenfunction of the Hamiltonian, 
the probability corresponding to the fcth box reads as 

E l*il 2 - (6) 

i£boxk 

One can introduce the f/tli moment of the box probability 
(frequently called generalized inverse participation ratio, 
GIPR), and its derivative: 

R i = E A 4 * S i = = E^fc ln ^' ( ? ) 

fc=i ^ fc=l 


The average of R q and S q follows a power-law behavior 
as a function of A = £, with exponent r q and a q : 


lim 

A->0 


In (Rq) 
In A 


= lim ^ . 

dq A->o {R q ) In A 


( 8 ) 


r q can be rewritten in the following form: 


T q = D q (q - 1) = d{q - 1) + A q , (9) 


where D q is the generalized fractal dimension, and A q is 
the anomalous scaling exponent. Employing a Legendre- 
transform on r g , we obtain the singularity spectrum, 

/(a) : 


f{a q ) = qot q - T q . (10) 

T q , a q , D q and A q are often referred to as multifractal 
exponents. 

According to recent resultP^ a symmetry relation ex¬ 
ists for a q and A q given in the form: 

A g = Ai_ g a q + ai- q = 2d (11) 


For numerical approaches one has to define the finite-size 
version of these MFEs at a particular value of disorder: 

( 12 ) 

D q ns (W, L, t) = * (13) 

where ens stands for ensemble averaging over the different 
disorder realizations. One may define typical averaged 
versions also: 





>J- 

' In A 
(In R q ) 
In A 


(14) 

(15) 


Similarly to a q and D q , A q or f q can be defined, which 
are called generalized multifractal exponents (GMFEs). 
Every GMFE approaches the value of the corresponding 
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MFE at the critical point, W = W Cl only in the limit A —> 
0. We would like to emphasize, that MFEs are defined 
through ensemble averaging in principle (see Eq. ©)> 
and ensemble and typical averaged MFEs are equal only 
in a range of q, q~ < q < defined by the two 

zeros of the singularity spectrum, f(a q _) = f(a q+ ) = 0. 
Therefore when in Sec. |IVB| we compute MFEs, we will 
use ensemble averaged quantities only. 

The choice of the investigated range of q is influenced 
by the following three effects. If q is large, the qth power 
in Eq. 0 enhances the numerical and statistical errors, 
leading to a noisy dataset. If q is negative with large 
absolute value, the relatively less precise small wave- 
function values dominate the sums in Eq. ([?]), which also 
results in a noisy dataset. These two effects together 
lead to a regime q m in < q < q m ax , where GMFEs be¬ 
have numerically the best. The third effect is coarse 
graining which suppresses the noise. For £ > 1 in an 
£ x £ x £ sized box positive and negative errors on the 
wave-functions can cancel each other. Moreover, in a 
box large and small wave-function amplitudes appear to¬ 
gether with high probability, and this way the relative 
error of a gk box probability is reduced. In other words 
coarse graining has a nice smoothing effect, which can 
help to widen the range of q that can be investigated. 

The renormalization flow of the AMIT has three fixed 
points: a metallic, an insulating and a critical one. In the 
metallic fixed point every state is extended with proba¬ 
bility one therefore the effective size of the states grows 
proportional to the volume, leading to D™ et = d. In 
the insulating fixed point every state is exponentially lo¬ 
calized, the effective size of a state does not change with 
changing system size, resulting in 7A™ S = 0 for q > 0, and 
D q ns = 00 for q < 0. Renormalization does not change 
the system at criticality, therefore it is scale independent, 
which means self-similarity. Therefore wave-functions are 
supposed to be multifractals, in other words generalized 
fractal^. 

Close to the critical point due to standard finite-size 
scaling arguments one can derive the following scaling 
laws for the exponents a q and D q defined above as: 


a q (W,L,£) =a q + ( q, 7 


D q (W,L,£) = D q + ] ^T q 


L 

L £ 


(16a) 

(16b) 


Equations (16aR(16b) can be summarized in one equa¬ 
tion: 


G q {W,L,£) = G q + ^G q 


L £ 

T'f 


(17) 


(L,£) on the left-hand side and |j on the right-hand 
side can be changed to ( L , A) and ( j, X\: 


G q (W,L,\) = G q + ^G q 



(18) 


Our central goal is to fit the above formulas to the nu¬ 
merically obtained data, where W c , v, y and G q appear 
among the fit parameters. This fit procedure will provide 
us the physically interesting quantities and their confi¬ 
dence intervals. In the next sections we present different 
methods for the finite-size scaling. 


A. finite-size scaling at fixed A 


At fixed A, G q in Eq. (18) can be considered as the 
constant term of Q q , therefore 


G q (W,L) 



(19) 


where the constant A has been dropped. Q q can be ex¬ 
panded with one relevant, g(w ), and one irrelevant oper¬ 
ator, r)(w), the following way by using w = W — W c : 

g q (eL*,nL-*) = g r q (pri) +r 1 L-yg i ; ( q l £) ( 20 ) 

All the disorder-dependent quantities in the above for¬ 
mula can be expanded in Taylor-series: 


e; 

*WI 

II 

( 21 ) 

s'; A) 

m r 

= 

2 — 0 

( 22 ) 

n e 

Urj 


g(w) = w + CiW 1 

r](w) = 1 + ^ diW 1 

(23) 


i-2 i =1 


The advantage of this method is, that in the Taylor-series 
only one variable appears, gL », therefore the number of 
parameters (including W Cl v and y) is n r + ni r + n p + n v + 
4, which grows linearly with the expansion orders. This 
method is very effective for computing W c , i/, and y, but 
since A is fixed, one cannot obtain the MFEs. In all cases 
we used A = 0.1, because it leads to excellent results in 
Ref. dH It seems, that it is small enough to capture the 
details of a wave-function, and it allows many different 
system sizes in the range of 20 < L < 100, which we 
investigated. This way we can also compare our results 
to those of Ref. M very well. 


B. Finite size scaling for varying A 

In order to take into account different values of A the 
scaling law given in Eq. © has to be considered. The 
expansion of g in © is 

g q ( e LKQ<’Kv'L- y, ,r 1 r y )=g r q (QLi,e£^ + 

+r I'L-yg'" (eLl,o&) +r 1 £~yg i ; . 
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According to Rodriguez et oiP-^ the most important ir¬ 
relevant term is the one containing the finite box size, t, 
therefore we took into account that one only. This leads 


to 


GqiW, L, i) = Gq + 

i^(s;(^,<A) + 

+ 

r,e-yg" (gLb,g^y (24) 

The Taylor expansions of the above functions are 

Gq {gL», 

n r i 

= VVu„A/.i/ (25) 

2=0 j—0 

Gq (gLKet") 

rii r i 

= ( 26 ) 

2=0 j — 0 

n e 

Tlrj 

q(w) = W + CiW % 

i=2 

rj(w) = 1 + ^ diW 1 (27) 
2=1 


The advantage of this method is, that it provides the 
MFE, G q , since it is one of the parameters to fit. There 
are many more data to fit compared to the fixed A case. 
Fixed A means that at a given system size one can use 
GMFEs obtained at a certain value of i - the one that 
leads to the desired A - , while in this case one can fit 
to GMFEs obtained at different values of t. However, 
these GMFEs are correlated, because they are the re¬ 
sults of the coarse graining of the same wave-functions 
with different sizes of boxes. During the fitting proce¬ 
dure one has to take into account these correlations, see 
Sec. |III C| Since the relevant and irrelevant scaling func¬ 
tions have two variables, gL* and g£", one has to fit 
a two-variable function with the number of parameters 
( TL r + 1) {jl r + 2)/2 + (fli r + 1) (Tlir T 2)/2 + Tip + 77^ + 3. We 
can see, that the number of parameters grows as ~ n^/ ir , 
instead of as ~ n r / ir as for fixed A. This makes the fitting 
procedure incorporating the correlations definitely much 
more difficult. 


and the correlation matrix of the numerically obtained 
data points by C, which can be computed numerically 
with a similar expression to the variance. With these 
notations y 2 reads as 

X 2 = - fi) (<?“% (Vi ~ /i)> (28) 

for more details see Ref. El If the data points are not 
correlated, C is a diagonal matrix, and the expression 
leads to the usual form: 

X 2 = E (29 ) 

i a i 

The number of degrees of freedom, Njf is the number of 
data points minus the number of fit parameters. A ratio 
X 2 /(A r rf/ — 1) « 1 means that the deviations from the best 
fit are of the order of the standard deviation (correlation 
matrix). The second criterion was that the fit has to be 
stable against changing the expansion orders, i.e. adding 
a few new expansion terms. From the fits that fulfilled 
the first two criteria we chose the simplest model, with 
the lowest expansion orders. Sometimes we also took into 
account the error bars, and we chose the model with the 
lowest error bar for the most important quantities (W c , v, 
etc...), if similar models fulfilled the first two criteria. 

The error bars of the best fit parameters were obtained 
by a Monte-Carlo simulation. The data points are results 
of averaging so due to central limit theorem, they have a 
Gaussian distribution. Therefore we generated Gaussian 
random numbers with parameters corresponding to the 
mean of the raw data points and standard deviation (or 
correlation matrix) of the mean, and then found the best 
fit. Repeating this procedure JVmc = 100 times provided 
the distribution of the fit parameters. We chose 95% 
confidence level to obtain the error bars. 


IV. RESULTS OF THE MFSS FOR THE 
ANDERSON MODELS IN THE WD SYMMETRY 
CLASSES 


C. General principles for the FSS fit procedures 


In this section we discuss the details of the methods 
and criteria we used during the MFSS. In order to fit 
the scaling law Eq. (19) and (24) we used the MINUIT 
library 22 *. To find the best fit to the data obtained nu¬ 
merically the order of expansion of Q T q lr , g and q must 
be decided by choosing the values of n r ,n,i r ,n e and n rl . 
Since the relevant operator is more important than the 
irrelevant one we always used n re i > rii r and n e > n v . To 
choose the order of the expansion we used basically three 
criteria. The first criterion we took into account was to 
check how close the ratio y 2 /(Ndf — 1) approached unity, 
where Ndf stands for the number of degrees of freedom. 
Let us denote the numerically obtained data points by y t . 
the fit function value at the ith parameter value by ft, 


With the numerical method described in Sec. ini we 
computed an eigenvector for every disorder realization 
of the Hamiltonian. From the eigenvectors every GMFE 
is computable, for the orthogonal and unitary class the 
I’lql 2 expression in Eq.([6| is trivial, and it means sum¬ 
mation for the spin-index for the symplectic class, since 
spatial behavior is in our interest. At fixed q exponents 
r q and A q are linear transforms of D q , so we used only 
the a q and D q GMFEs for the MFSS. We investigate the 
range — 1 < q < 2, because GMFEs behave the best in 
this regime for the reasons described in Sec. |III| 


A. Results of the MFSS at fixed A = 0.1 

The typical behavior of the GMFEs is presented in 
Fig-B In all cases there is a clear sign of phase transition: 
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FIG. 1: Dots are the raw data for different GMFEs in the conventional WD symmetry classes. Red line is the best fit obtained 
by MFSS. Insets are scaling functions on a log-log scale, after the irrelevant term was subtracted. Error bars are shown only 
on the large figures, in order not to overcomplicate the insets. 


With increasing system size the GMFEs tend to oppo¬ 
site direction on both sides of their crossing point. Note 
that there is no well-defined crossing point due to the 


irrelevant term in Eq. (201. Applying the MFSS method 
described in Sec. Ill A with the principles of Sec. IIIC 
to the raw data leads to a well- fitting function, see red 
lines in Fig. [I] After the subtraction of the irrelevant 
part from the raw data, plotting it as a function of gL» 
results a scaling-function also, see insets of Fig. [T] 

The MFSS provided us the critical point, W c , the crit¬ 
ical exponent, v, and the irrelevant exponent, y at every 
investigated values of q, the results are given in Fig. [2] 
The parameters of the critical point correspond to the 
system itself, therefore it should not depend on the quan¬ 
tity we used to find it. In other words, it should be in¬ 
dependent of q, the averaging method and the GMFE 
we used. From Fig. [2] it is clear that this requirement 
is fulfilled very nicely. There is a small deviation for the 
irrelevant exponent, y , obtained from a typ at q = —1 and 
q = —0.75 in the unitary and symplectic class, but since 
y describes the subleading part, it is very hard to deter¬ 
mine, and we cannot exclude some sort of underestimati- 
ion of the error bar of this exponent. Another interesting 
feature of the results is that the error bars get larger as 
q goes above 1. As written in Sec. |III C] large q enhances 
the errors through the gth power in Eq. leading to 
bigger error bars. A similar effect can be seen around 
q ss —1, where the relatively less precise small wave- 
function values dominate the sums in Eq. ([T]), which can 
also contribute to the deviation of y obtained from a typ 
in this regime. These two effects together lead to our 
investigated interval — 1 < q < 2, where GMFEs behave 
the best. The results are strongly correlated, since they 
were obtained from the same wave-functions, therefore 
they cannot be averaged. We chose a typical g-point for 
every symmetry class to describe the values of the critical 
parameters, see Tab. [TT] 

In the orthogonal class the critical parameters are 
in excellent agreement with the most recent high 
precision results of Rodriguez et al P^, W^ od = 
16.517 (16.498..16.533), v% x d = 1.612 (1.593..1.631) and 


VrocI = 1-67 (1.53..1.80), obtained from d 0 with the same 
method (fixed A). This agreement verifies our numerics 
and fit method, and makes it reliable for the other two 
universality classes. 

In the unitary class the critical parameters match 
with the results of Slevin and Ohtsukf^, Sle = 

18.375 (18.358..18.392) and v u sie = 1.43 (1.37..1.49), ob¬ 
tained by transfer matrix method (they did not published 
the value of the irrelevant exponent). They used mag¬ 
netic flux <!> = j , while we used $ = |, and according 
to Drose et alW- 1, depends on the applied magnetic 

flux. However, in Fig.2. of Ref. [23] it can be seen that 
the critical points at $ = \ and <I> = I are very close 
to each other, hence the agreement between our critical 
point and the result of Slevin and Ohtsuki. 

In the symplectic class the critical parameters agree 
more or less with the results of Asada et a/P^O, W r f Aso = 
20.001 (19.984..20.018), u% sa = 1.375 (1.359..1.391) 
and y Asa = 2.5 (1.7..3.3), obtained by transfer matrix 
method. However, the difference does not seem to be very 
large, our critical point is considerably different, even 
though we used exactly the same model. Due to big¬ 
ger computational resources we could investigate much 
bigger system sizes than they did, therefore it is possi¬ 
ble that they underestimated the role of the irrelevant 
scaling, resulting in a somewhat higher critical point. 

The critical points are higher in the unitary and in the 
symplectic class, than in the orthogonal class, showing 
that broken time-reversal or spin-rotational symmetry re¬ 
quires more disorder to localize wave-functions. Since the 
value of the critical point in the unitary and symplectic 
class can be influenced by the strength of the applied 
magnetic flux and spin-orbit coupling, the relationship 
between W^ x and W^ x probably depends on these two 
parameters. However, because of their close value of the 
critical exponents, and v sx are the same within our 
confidence interval, and the following relation appears: 
v ox > v ux > v sx . The situation for the irrelevant ex¬ 
ponent is similar namely, that they are the same within 
error bar, but y ox seems to be slightly higher than y ux , 
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FIG. 2: Critical parameters of the Anderson models in WD classes obtained by MFS at fixed A = 0.1. First row corresponds 
to the orthogonal class, second row corresponds to the unitary class, and third row corresponds to the symplectic class. 


class 

exp 

w c A 


y A 

N df 

X 

TlrTlirTl gTlri 

ort 

~ ens 

O0.6 

16.524 (16.511..16.538) 

1.598 (1.576..1.616) 

1.763 (1.679..1.842) 

172 

176 

3 2 10 

uni 

~ ens/typ 

a 0 

18.373 (18.358..18.386) 

1.424 (1.407..1.436) 

1.633 (1.516..1.751) 

198 

179 

4 2 10 

sym 

f)typ 

-^-0.25 

19.838 (19.812..19.869) 

1.369 (1.305..1.430) 

1.508 (1.309..1.743) 

171 

151 

4 2 10 


TABLE II: Result of the MFSS at fixed A = 0.1 for the selected values of q. 


which is a bit higher than y sx . 


B. Results of the MFSS at varying A 


As mentioned in Sec. Ill C GMFEs obtained by typical 
averaging are equal to ensemble -averaged GMFEs only in 
a range of q, q~ < q < q+. Since we intend to compute 
the MFEs also, we restrict our analysis to ensemble aver¬ 
aged GMFEs, and drop the label ens from the notation. 


We fit the formula Eq. (24 (to the raw data. To do 


that, we choose a range of box size £, which is used for 
the MFSS. We always use the widest range of £, that re¬ 
sults in convergence, X 2 /{^df — 1) ~ 1. We find that for 
our dataset for different values of q for a q or D q different 
ranges of l were the best. We used minimal box sizes 
i rn i n = 2 or £ m i n = 3 and maximal box sizes correspond¬ 
ing to \ max = 0.1 or X max = 0.066. At a 0 A and a 0 . 6 
the fitting method sometimes suffered from convergence 
troubles and resulted in large error bars, because these 
points are close to the special case of q = 0.5 where, by 
definition, op.5 = d. Artifacts from this regime were also 



















































































reported in Ref. HU so we decided not to take into ac¬ 
count these points for a. We tried several combinations of 
£min, A max and expansion orders in the symplectic class 
for ai .75 and a^, but none of them resulted in stable 
fit parameters. Therefore values computed from these 
points are also missing from our final results, which are 
visible in Fig. [3] The results are independent of q and 
the GMFE we used, similar to the fixed A method. In 
Sec. IV A we already saw that according to the arguments 
of Sec. Ill C| error bars get bigger, if q grows beyond 1. 
This phenomenon is more amplified here, especially for 
values coming from fits for a q , but larger error bars on 
values corresponding to D q are present on a moderate 
level also. Since Fig .8 of Ref. HU shows results for this 
regime only for values corresponding to A q , which is a 
linear transform of D qi we can compare their results only 
to ours corresponding to D q . One can see that our error 
bars are similar, even though there are differences prob¬ 
ably due to the fact that they used system sizes up to 
L = 120, which was not possible for us, mainly because 
of the long runtime and large memory usage for the sym¬ 
plectic model. They also use £ m i n = 1 and £ m ; n = 2, 
while imin = 1 was never suitable for our dataset. We do 
not know the precise origin of this behavior, but we have 
a few possible explanations. We experience that larger 
system sizes allow a wider range of £ to be used. We 
have smaller system sizes than Ref. HU and fewer sam¬ 
ples for the largest systems sizes. Noise also gets bigger 
as £ decreases, because of the smoothing effect of box¬ 
ing described in Sec. |IIIC[ which can also explain partly 
our experience. Another important difference is that in 
Eq. (37) of Ref. HU the authors use an expression in the 
expansion of the scaling function, which is proportional 
to the square of the irrelevant term, (rj £~ v ) 2 . According 
to our experience the inclusion of this term produced no 
improvement in the scaling analysis, so we use the scaling 
function described in Eq. [24] Such a difference might be 
explained again by our different dataset. 


C. Analysis of the multifractal exponents 


MFSS for varying A provided us the MFEs in all WD 
classes, which are listed in Tab.[Vj and depicted in Fig. [4] 
For the orthogonal class one can find matching results 
with the listed MFE-s in Ref. HU Since the precise values 
of the MFEs in three dimensions were determined first in 
Ref. HU for the orthogonal class only, the lack of reliable 
analytical and numerical results for the other symmetry 
classes makes our results more important. The most con¬ 
spicuous thing in Fig. [4] is that curves for different sym¬ 
metry classes are very close to each other, they are almost 
indistinguishable at the first sight. This shows that the 
broken time-reversal or spin rotational symmetry has a 
very small effect on the MFEs in three dimensions. Tak¬ 
ing a closer look (or from Tab. [V} one can see that the 
curve of D q and a q are the steepest in the symplectic, 
the second steepest in the unitary, and the less steep in 
the orthogonal class. From Tab. [V] it is also clear that 
at most of the q values there is a significant difference 
between the MFEs of different symmetry classes. 


There are no critical states in the two dimensional or¬ 
thogonal clasf^, but one can find values of ao for the 
two dimensional unitary class (Integer Quantum Hall), 
a 02 D = 2.2596 ± O.OOOl^, and symplectic class, aQ 2D = 
2.172 ± O.OOSP^. Comparing the difference between these 
exponents in two dimensions we get a^ 2D — a 02 d = 
0.0876 ± 0.0024, while our result for three dimensions is 
a 0 3D ~ a 0 3D = —0.03 ±0.015. There is about a factor of 
3 between the magnitude of these values, and even their 
sign is opposite, which shows very different effect of pres¬ 
ence or absence of spin rotational symmetry in different 
dimensions. 


As written in Sec. |IV A[ the results for different values 
of q are strongly correlated, therefore we chose one of 
them with the lowest error bars that represents well the 
results for that universality class. 


The critical parameters listed in Tab. m are in a very 
nice agreement with our previous results for the fixed 
method of A = 0.1, see Sec. |IV A[ and also with the 
results of Refs. HUHU Comparing the critical param¬ 
eters for the orthogonal case with the results of Ro¬ 
driguez et obtained by the same method, W° Rod = 
16.530 (16.524..16.536), v% od = 1.590 (1.579..1.602), we 
see a nice agreement again. Moreover these results are 
more accurate with this method compared to the fixed A 
method, leading to (for y° and y u only almost) signif¬ 
icantly different critical exponents and irrelevant expo¬ 
nents for the different WD classes, v° > v u > v s and 

y° >y u > y s - 


We tested the symmetry relation Eq. ( pTTj ) for a q and 
A g , the results are listed in Tab. [V] and depicted in 
Fig-0 The symmetry relation is fulfilled in the range 
—0.25 < q < 1.25 (in the symplectic class only for 
—0.25 < q < 1), and small deviations are visible out¬ 
side this interval. In this regime error bars are growing 
very large, coming mainly from the large errors of a q >i .5 
and Dq>i 5 . Similar effects were already seen for the crit¬ 
ical parameters in Fig. [3] It is really hard to estimate the 
correct error bars in this large q case, and the deviations 
from symmetry are small, therefore we believe that dif¬ 
ferences appear only because of slightly underestimated 
error bars of a q > 1.5 and D q >i 5 . All in all we find nu¬ 
merical results basically matching with Eq. (11). 


Assuming, that A q is an analytic function of q , and 
using the symmetry relation, Eq. (11), one can expand 
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FIG. 3: Critical parameters of the Anderson models in WD classes obtained by two-variable MFFS with varying A. First row 
corresponds to the orthogonal class, second row corresponds to the unitary class, and third row corresponds to the symplectic 
class. 


class 

exp 

W c 

V 

y 

N df 


Th 1- Tl i t Th q Tt tj 

ort 

do 

16.524 (16.513..16.534) 

1.595 (1.582..1.609) 

1.749 (1.697..1.786) 

241 

267 

3 2 10 

uni 

D0.1 

18.371 (18.363..18.380) 

1.437 (1.426..1.448) 

1.651 (1.601..1.707) 

275 

232 

4 2 10 

sym 

do 

19.836 (19.831..19.841) 

1.383 (1.359..1.412) 

1.577 (1.559.. 1.595) 

361 

352 

3 2 10 


TABLE III: Critical parameters of the Anderson models in the WD symmetry classes obtained by two-variable MFSS with 
varying A. 


A q in Taylor series around q = \: 


\ = £ 


Ck q 


k—0 


2 k oo 


k=0 


= E Cfc 9(5-!)+ 7 = 


k =0 *=0 w 

OO 

(g(l-<?))\ 


k—i 


(30) 


k=1 


where the condition Ao = Ai = 0 enforced by the defini¬ 
tion of A q (see Eq. ([9])) was used in the last step, leading 
to k = 1 as the lower bound for the summation. Similar 
expression can be derived for a q by using the connection 
a q = d + A q derived from Eqs. |8j)-(|9j): 

OO 

a q = d+ (1 - 2g)^a fc (g(l - q)) k ~ r , (31) 

k=1 

where a& = kdk , and a\ = d\ = ckq — d. One can 
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FIG. 4: MFEs of the Anderson models in the WD universality 
classes. Corresponding data are listed in Tab. [V] 
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FIG. 5: Test for symmetry relation Eq. (Ill in the WD sym¬ 
metry classes. Points are shifted horizontally a little bit for 
better visualization. Only the range q > 0.5 is visible because 
expression a q -\-ai- q (A 9 —Ai_ g ) is symmetric (antisymmet¬ 
ric) for q = 0.5. 


obtain the dk and a*, coefficients by fitting the expres¬ 
sions Eq. (30) (31). We used only the range q < 1.25, 
because beyond this regime error bars are growing ex¬ 
tremely large, and the re a re small deviations from the 
symmetry relation Eq.|i~tj) also. We plotted and 

°i- 2 q Fig- [®] make the presence of higher-order 
terms of the expansion visible. 

We fit expressions Eq. (30) (31) up to third order in 
all cases, the resulting expansion coefficients are listed in 
Tab. [TV] From the data listed one can see that the ex¬ 
pansion coefficients fulfill the relation ak = kdk■ How¬ 
ever ct q and A q were obtained from the same wave- 
functions, they are results of completely independent fit- 
procedures. Therefore the fact, that they satisfy the 
equation a*, = kdk further confirms our result for their 
value listed in Tab. |Vl for q < 1.25 and shows the consis¬ 
tency of the MFSS. 

As one would expect for expansion coefficients, dk and 
Ofc show decreasing behavior as k grows. Only d\ and 


1.25 
1.2 
1.15 
1.1 

Tiff 1.05 

iU 1 

0.95 
0.9 
0.85 
0.8 

-2 -1.5 -1 -0.5 0 0.5 


orthogonal 

unitary 

symplectic 

i V ■ 

3=^ 

i V ■ 

5 


1.2 

1.15 

orthogonal i— ■—i 

unitary ■—♦—< _ . 

1.1 

symplectic >—-—' jp-'' 

1.05 


1 

■ ' 

0.95 



-T" 


-2 -1.5 -1 -0.5 0 0.5 


9(1 - 9) 


9(1 -9) 


FIG. 6: Dots and error bars are numerical values for the cor¬ 
responding quantities, and , for the WD symme¬ 

try classes. Lines are the best fits. Several points are shifted 
horizontally a bit for better viewing. 



ort 

uni 

sym 

di 

1.044 (1.041..1.047) 

1.097 (1.095..1.098) 

1.123 (1.122..1.125) 

d 2 

0.095 (0.085..0.105) 

0.096 (0.091..0.100) 

0.088 (0.084..0.093) 

d 3 

0.018 (0.011..0.025) 

0.017 (0.014..0.020) 

0.014 (0.010..0.017) 

ai 

1.045 (1.042..1.048) 

1.099 (1.096.. 1.102) 

1.124 (1.123..1.126) 

a 2 

0.182 (0.168..0.195) 

0.185 (0.174..0.197) 

0.185 (0.179..0.191) 

a 3 

0.044 (0.035..0.053) 

0.043 (0.035..0.050) 

0.044 (0.038..0.049) 


TABLE IV: Expansion coefficients of Eqs. ( |30| )-( |3Tj ) obtained 
by a fit depicted in Fig. [6] 


(ii are significantly different for the different symmetry 
classes, while d 2 , cfe, et 2 and 03 are the same within er¬ 
ror bars. Their real value is probably different, but the 
relative error of the expansion coefficients naturally in¬ 
creases as k grows, leading to indistinguishable values for 
the different symmetry classes for k > 2 . 

Wegner computed analytically 2 -^ the value of A q with 
e expansion using nonlinear <r-model up to fourth-loop 
order for the orthogonal and the unitary symmetry class, 
resulting an expansion in dimensions d = 2 + e for £ -C 
lEE % 


A° q = g(l - q)e + ^ q(q - l)(g 2 - q + 1) £ 4 + 0(e 5 ) 

= (e~ ^ 4 ) 9(1 - 9 ) + ^£ 4 ( 9(1 - 9)) 2 + 0(e 5 ) 

(32) 

= y|«(l - 9 ) - §C( 3 ) £ 2 ( 9(1 - q)) 2 + O(si) (33) 

Even though £< 1 should hold, one can try to extrapo¬ 
late to three-dimensions by inserting e = 1. This leads to 
d? « 0.699, d§ « 0.301, df « 0.707 and d^ « -0.451. 
As one can see, these values are rather far from our 
numerical results, but this is not surprising for an £- 
expansion at e = 1. These results capture well the ten¬ 
dency at least that d° is slightly smaller, than dd[. On the 
other hand it leads to d 2 and d 2 having opposite sign, 
which is highly inconsistent with our numerical results. 
It is interesting that the first-loop term, which is propor- 
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tional to e and leads to parabolic A q , results in = 1 
and a® = ao — d = 1, which are very close to our numer¬ 
ically measured values. In this sense parabolic approx¬ 
imation is better for the orthogonal class, as compared 
to the fourth-loop order approximation. If higher-order 
terms were obtained, or if A q were expanded by using an¬ 
other approach, our coefficients could provide relatively 
accurate values as compared with analytical results. 


V. SUMMARY 


In this paper we examined the three-dimensional An¬ 
derson models belonging to the conventional WD symme¬ 
try classes with the help of multifractal finite-size scal¬ 
ing using two methods: a simpler method for fixed A 
leading to a single-variable scaling function, and a more 
complicated one for varying A resulting in a two-variable 
scaling function. Both methods confirmed the presence 
of multifractality in all three symmetry classes, and we 
obtained critical parameters listed in Tabs. [TT] and III 
in aggreement with each other and with previous results 
known from the literature. The more complicated vary¬ 
ing A method provided more precise values for the critical 


parameters, listed in Tab. Ill and significantly different 
critical exponents for the different WD symmetry classes. 

Applying the method of varying A we also calculated 
the multifractal exponents, that basically fulfill the ex¬ 
pected symmetry relation Eq. , small deviations were 


detected for large q -values probably due to slightly under¬ 
estimated error bars, in Fig.[4]one can see that the MFEs 
of different symmetry classes are very close to each other, 
but Fig.[6]or Tab. |V| shows significant differences between 
them for most of the values of q. We compared the differ¬ 
ence of ao in the unitary and symplectic class to available 
results in two dimensions, and we found completely dif¬ 
ferent relation between the two and three dimensional 
cases. We expanded the MFEs in terms of the variable 
q{ 1 — q ), and determined the expansion coefficients up 
to third order numerically. The expansion coefficients 


of Eq. (30) (31) fulfill the expected relation a*, = kdk 
giving a further confirmation for the validity of our re¬ 
sults for the MFEs listed in Tab. [V] We also compared 
the numerical results to available analytical estimates, 
and found in some cases similar, but in other cases op¬ 
posite qualitative behavior for expansion coefficients for 
the orthogonal and the unitary classes. Nevertheless, we 
believe that the numerical precision of our results should 
be used as tests for future renormalization or other type 
of expansion approximations. Therefore our results await 
analytical comparison. 
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q 

class 

a Q 

D q 

/w 


y q +a 1 


^<2 

-Ai- 



ort 

5.555 

(5.490. 

5.626) 

3.926 

(3.914. 

.3.938) 

2.297 

(2.338. 

2.250) 

6.275 

(6.042. 

6.661) 

-0.102 ( 

-0.218 

.0.000) 

-i 

uni 

5.671 

(5.629. 

5.707) 

3.970 

(3.966. 

.3.976) 

2.269 

(2.303. 

2.245) 

6.331 

(6.215. 

6.444) 

-0.130 (- 

0.195. 

-0.062) 


sym 

5.751 

(5.690. 

5.799) 

4.001 

(3.994. 

.4.010) 

2.251 

(2.298. 

2.222) 

6.379 

(6.197. 

6.584) 

-0.134 (- 

0.237. 
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ort 
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(5.187. 

5.267) 
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(3.708. 

.3.722) 
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2.564) 
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(5.988. 

6.353) 

-0.035 ( 

-0.094 

.0.032) 

-0.75 
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(5.317. 

5.349) 
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(3.749. 

.3.754) 
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(2.573. 
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6.176 

(6.131. 

6.239) 

-0.062 (- 
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-0.025) 


sym 

5.406 

(5.387. 

5.430) 

3.773 

(3.770. 

.3.777) 

2.549 

(2.558. 

2.537) 

6.221 

(6.113. 

6.349) 

-0.060 (- 

0.114. 

-0.023) 


ort 
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(4.856. 

4.896) 
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(3.488. 

.3.496) 

2.800 

(2.803. 

2.796) 
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(5.959. 

6.149) 

-0.008 ( 

-0.045 

.0.025) 

-0.5 

uni 

4.975 

(4.958. 

4.994) 

3.517 

(3.512. 

.3.521) 

2.788 

(2.789. 

2.785) 

6.103 

(6.000. 

6.167) 

-0.009 ( 

-0.025 
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sym 
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(5.019. 

5.039) 

3.532 

(3.531. 

.3.534) 

2.784 

(2.787. 

2.781) 
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(6.039. 

6.206) 

-0.019 (- 

0.041. 

-0.001) 


ort 

4.488 

(4.477. 

4.499) 
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.3.255) 
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sym 
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ort 
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-0.005 

.0.003) 

0.2 

uni 

3.678 

(3.673. 

3.684) 

2.778 

(2.777. 

.2.780) 

2.958 

(2.956. 

2.961) 

5.999 

(5.987. 

6.011) 

-0.000 ( 

-0.004 

.0.004) 


sym 

3.693 

(3.691. 

3.695) 

2.772 

(2.772. 

.2.773) 

2.956 

(2.955. 

2.957) 

5.999 

(5.995. 

6.004) 

0.001 (- 

0.000. 

0.003) 


ort 

3.541 

(3.534. 

3.547) 

2.734 

(2.733. 

.2.737) 

2.936 

(2.933. 

2.939) 

6.000 

(5.987. 

6.012) 

-0.001 ( 

-0.006 

.0.003) 

0.25 

uni 

3.569 

(3.563. 

3.575) 

2.721 

(2.720. 

.2.722) 

2.933 

(2.931. 

2.935) 

5.999 

(5.987. 

6.011) 

-0.000 ( 

-0.002 

.0.001) 


sym 

3.579 

(3.577. 

3.581) 

2.715 

(2.714. 

.2.715) 

2.931 

(2.930. 

2.932) 

5.997 

(5.992. 

6.001) 

0.001 (- 

0.001. 

0.003) 


ort 

3.436 

(3.430. 

3.441) 

2.681 

(2.678. 

.2.684) 

2.907 

(2.903. 

2.911) 

6.001 

(5.991. 

6.012) 

-0.001 ( 

-0.006 

.0.004) 

0.3 

uni 

3.459 

(3.453. 

3.464) 

2.665 

(2.664. 

.2.666) 

2.903 

(2.900. 

2.905) 

5.999 

(5.987. 

6.010) 

-0.000 ( 

-0.002 

.0.001) 


sym 

3.465 

(3.462. 

3.467) 

2.657 

(2.656. 

.2.658) 

2.899 

(2.898. 

2.901) 

5.995 

(5.991. 

6.000) 

0.001 (- 

0.001. 

0.003) 


ort 


- 


2.573 

(2.570. 

.2.577) 


- 



- 


-0.001 ( 

-0.006 

.0.004) 

0.4 

uni 


- 


2.551 

(2.550. 

.2.553) 


- 



- 


-0.000 ( 

-0.002 

.0.002) 


sym 


- 


2.542 

(2.540. 

.2.543) 


- 



- 


0.001 (- 

0.001. 

0.003) 


ort 


3 (3..3) 

2.466 

(2.459. 

.2.471) 

2.733 

(2.730. 

2.736) 


6 (6..6) 

0 

(0..0) 


0.5 

uni 


3 (3..3) 

2.439 

(2.437. 

.2.441) 

2.719 

(2.719. 

2.721) 


6 (6..6) 

0 

(0..0) 



sym 


3 (3..3) 

2.427 

(2.425. 

.2.429) 

2.714 

(2.712. 

2.715) 


6 (6..6) 

0 

(0..0) 



ort 


- 


2.358 

(2.352. 

.2.366) 


- 



- 


0.001 (- 

0.004. 

0.006) 

0.6 

uni 


- 


2.327 

(2.325. 

.2.329) 


- 



— 


0.000 (- 

0.002. 

0.002) 


sym 


- 


2.314 

(2.311. 

.2.317) 


- 



- 


-0.001 ( 

-0.003 

.0.001) 


ort 

2.566 

(2.561. 

2.571) 

2.252 

(2.242. 

.2.263) 

2.472 

(2.466. 

2.479) 

6.001 

(5.991. 

6.012) 

0.001 (- 

0.004. 

0.006) 

0.7 

uni 

2.540 

(2.535. 

2.545) 

2.217 

(2.214. 

.2.220) 

2.443 

(2.438. 

2.448) 

5.999 

(5.987. 

6.010) 

0.000 (- 

0.001. 

0.002) 


sym 

2.530 

(2.528. 

2.532) 

2.203 

(2.199. 

.2.207) 

2.432 

(2.429. 

2.435) 

5.995 

(5.991. 

6.000) 

-0.001 ( 

-0.003 

.0.001) 


ort 

2.459 

(2.454. 

2.465) 

2.198 

(2.186. 

.2.209) 

2.394 

(2.387. 

2.401) 

6.000 

(5.987. 

6.012) 

0.001 (- 

0.003. 

0.006) 

0.75 

uni 

2.430 

(2.424. 

2.436) 

2.163 

(2.159. 

.2.168) 

2.363 

(2.358. 

2.369) 

5.999 

(5.987. 

6.011) 

0.000 (- 

0.001. 

0.002) 


sym 

2.417 

(2.415. 

2.419) 

2.148 

(2.143. 

.2.156) 

2.350 

(2.347. 

2.353) 

5.997 

(5.992. 

6.001) 

-0.001 ( 

-0.003 

.0.001) 


ort 

2.354 

(2.347. 

2.360) 

2.147 

(2.135. 

.2.157) 

2.312 

(2.304. 

2.319) 

5.998 

(5.985. 

6.011) 

0.001 (- 

0.003. 

0.005) 

0.8 

uni 

2.320 

(2.314. 

2.326) 

2.111 

(2.099. 

.2.125) 

2.278 

(2.271. 

2.286) 

5.999 

(5.987. 

6.011) 

0.000 (- 

0.004. 

0.004) 


sym 

2.307 

(2.304. 

2.309) 

2.095 

(2.090. 

.2.100) 

2.264 

(2.261. 

2.267) 

5.999 

(5.995. 

6.004) 

-0.001 ( 

-0.003 

.0.000) 


ort 

2.146 

(2.135. 

2.159) 

2.046 

(2.029. 

.2.060) 

2.136 

(2.124. 

2.149) 

5.995 

(5.978. 

6.014) 

0.001 (- 

0.002. 

0.003) 

0.9 

uni 

2.107 

(2.097. 

2.117) 

2.009 

(1.991. 

.2.025) 

2.097 

(2.087. 

2.108) 

5.997 

(5.981. 

6.014) 

-0.000 ( 

-0.002 

.0.002) 


sym 

2.092 

(2.088. 

2.099) 

1.988 

(1.981. 

.1.997) 

2.082 

(2.077. 

2.088) 

6.005 

(5.998. 

6.013) 

-0.001 ( 

-0.002 

.0.000) 


ort 

1.948 

(1.930. 

1.963) 


a l 



«l 


5.991 

(5.965. 

6.012) 

0 

(0..0) 


1 

uni 

1.905 

(1.886. 

1.925) 


a l 



a l 


6.000 

(5.974. 

6.026) 

0 

(0..0) 



sym 

1.886 

(1.877. 

1.896) 


a 1 



ai 


6.010 

(5.999. 

6.023) 

0 

(0..0) 



ort 

1.520 

(1.508. 

1.535) 

1.727 

(1.715. 

.1.738) 

1.477 

(1.418. 

1.551) 

6.009 

(5.985. 

6.034) 

-0.001 ( 

-0.006 

.0.003) 

1.25 

uni 

1.473 

(1.442. 

1.499) 

1.688 

(1.660. 

.1.708) 

1.422 

(1.391. 

1.457) 

6.036 

(5.995. 

6.073) 

0.006 (- 

0.002. 

0.011) 


sym 

1.437 

(1.424. 

1.450) 

1.644 

(1.634. 

.1.655) 

1.371 

(1.338. 

1.409) 

6.044 

(6.027. 

6.061) 

0.004 (0.001..0.007) 


ort 

1.185 

(1.161. 

1.206) 

1.534 

(1.518. 

.1.550) 

1.007 

(0.912. 

1.079) 

6.061 

(6.017. 

6.102) 

0.005 (- 

0.009. 

0.019) 

1.5 

uni 

1.096 

(1.073. 

1.124) 

1.468 

(1.453. 

.1.483) 

0.958 

(0.836. 

1.017) 

6.072 

(6.031. 

6.118) 

0.009 (- 

0.004. 

0.025) 


sym 

1.060 

(1.044. 

1.080) 

1.450 

(1.437. 

.1.465) 

0.889 

(0.827. 

1.011) 

6.090 

(6.063. 

6.118) 

0.024 (0.015..0.034) 


ort 

0.920 

(0.889. 

0.949) 

1.372 

(1.349. 

.1.395) 

0.590 

(0.422. 

0.818) 

6.145 

(6.076. 

6.216) 

0.029 (- 

0.001. 

0.058) 

1.75 

uni 

0.841 

(0.814. 

0.873) 

1.301 

(1.273. 

.1.329) 

0.479 

(0.459. 

0.529) 

6.175 

(6.130. 

6.222) 

0.041 (0.015..0.065) 


sym 

no stability 

1.262 

(1.242. 

.1.290) 







0.050 (0.030..0.077) 


ort 

0.719 

(0.683. 

0.754) 

1.231 

(1.203. 

.1.256) 

0.190 

-0.068 

..0.727) 

6.274 

(6.173. 

6.380) 

0.083 (0.031..0.132) 

2 

uni 

0.622 

(0.583. 

0.690) 

1.173 

(1.147. 

.1.205) 

0.131 

(0.039. 

0.230) 

6.293 

(6.212. 

6.396) 

0.113 (0.076..0.154) 


sym 

no stability 

1.118 

(1.099. 

.1.167) 







0.120 (0.083..0.184) 


TABLE V: MFE a q , D q and f(a q ), and values for the corresponding symmetry relation Eq. (Ill obtained for the WD symmetry 
classes. 































